c->cd(1);
h1->ProfileX("_pfx1", 0, 400);
_pfx1->Draw("same");
TF1 *rt = new TF1("rt", "pol5");
// rt->SetParLimits(0, -0.2, 0.3);
// rt->FixParameter(5, 0.);
_pfx1->Fit(rt, "QBR", "", -10., 110.);
for(int i = 0; i <= 5; i++){ cout << rt->GetParameter(i) << (i<5?", ":"\n"); }
cout << "r(0) = " << rt->Eval(0.) << " r(66) = " << rt->Eval(66.) << endl;

c->cd(5);
h5->ProfileX("_pfx5", 0, 400);
_pfx5->Draw("same");
TF1 *rt5 = new TF1("rt5", "pol5");
// rt5->SetParLimits(0, -0.2, 0.3);
// rt->FixParameter(5, 0.);
_pfx5->Fit(rt5, "QBR", "", -10., 110.);
for(int i = 0; i <= 5; i++){ cout << rt5->GetParameter(i) << (i<5?", ":"\n"); }
cout << "r(0) = " << rt5->Eval(0.) << " r(66) = " << rt5->Eval(66.) << endl;

c->cd(9);
h9->ProfileX("_pfx9", 0, 400);
_pfx9->Draw("same");
TF1 *rt9 = new TF1("rt9", "pol3");
// rt9->SetParLimits(0, -0.2, 0.9);
_pfx9->Fit(rt9, "QBR", "", -30., 300.);
for(int i = 0; i <= 3; i++){ cout << rt9->GetParameter(i) << (i<3?", ":"\n"); }
cout << "r(0) = " << rt9->Eval(0.) << " r(180) = " << rt9->Eval(180.) << endl;

c->cd(13);
h13->ProfileX("_pfx13", 0, 400);
_pfx13->Draw("same");
TF1 *rt13 = new TF1("rt13", "pol3");
// rt13->SetParLimits(0, -0.2, 0.9);
_pfx13->Fit(rt13, "QBR", "", -30., 300.);
for(int i = 0; i <= 3; i++){ cout << rt13->GetParameter(i) << (i<3?", ":"\n"); }
cout << "r(0) = " << rt13->Eval(0.) << " r(180) = " << rt13->Eval(180.) << endl;

c->cd(14);
h14->ProfileX("_pfx14", 0, 400);
_pfx14->Draw("same");
TF1 *rt14 = new TF1("rt14", "pol3");
// rt14->SetParLimits(0, -0.2, 0.9);
_pfx14->Fit(rt14, "QBR", "", -30., 300.);
for(int i = 0; i <= 3; i++){ cout << rt14->GetParameter(i) << (i<3?", ":"\n"); }
